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Background: Uncertainty in the operating temperatures of Oklo reactor zones impacts the precision of bounds derived for 
time variation of the fine structure constant a. Improved 176 Lu/ 175 Lu thermometry has been discussed but its usefulness 
may be complicated by photo excitation of the isomeric state 176m Lu by 176 Lu(7, 7') fluorescence. 

Purpose: We calculate prompt, delayed and equilibrium 7-ray fluxes due to fission of 235 U in pulsed mode operation of Oklo 
zone RZ10. 

Methods: We use Monte Carlo modeling to calculate the prompt flux. We use improved data libraries to estimate delayed 
and equilibrium spectra and fluxes. 

Results: We find 7-ray fluxes as a function of energy and derive values for the coefficients A 7i y that describe burn-up of 176 Lu 
s ! ' through the isomeric 176m Lu state. 

Conclusion: The contribution of the (7, 7') channel to the 176 Lu/ 175 Lu isotopic ratio is negligible in comparison to the neutron 
^ 1 burn-up channels. Lutetium thermometry is fully applicable to analyses of Oklo reactor data. 

PACS numbers: 06.20.Jr, 07.05.Tp, 25.20.Dc, 25.85.Ec, 28.20.Gd, 28.41. Kw 
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1. Introduction. 

— — Studies of the 235 U fission product isotopic ratios from Oklo jjO] have been undertaken by many groups investigating 
whether the fine structure constant a has changed over the 2 GY period since the reactors operated. As first pointed 
out by Shlyachter the samarium isotopic ratios are sensitive to the value of a through the overlap of the 149 Sm 
Eo=97.3-meV neutron resonance with the thermal and epithermal portions of the neutron flux in the reactor. While 
^vq , the majority of Oklo analyses [3j-|7j], have been consistent with no shift in the resonance energy, and therefore no 
t-H ■ change in a, a change has been argued for from astronomical observations Q. 

All Oklo analyses make assumptions about the operating temperatures of the reactors. But there is as yet no 
■ agreement on what these temperatures actually were. Utilizing the 176 Lu/ 175 Lu isotope ratio method to determine 
7— I I temperatures was recently revisited by Gould and Sharapov 9] . The method is based on the temperature dependence 
of the large thermal neutron capture cross section of 176 Lu (natural abundance 2.599 % and on knowing with 

certainty the (small) ground state branching ratio for thermal neutron capture on the more abundant lutetium isotope, 
175 Lu (natural abundance 97.401 %). The dominant capture branch a^} 5 leads to a short lived isomeric state in 176m Lu, 
while only a minor branch <rf 75 leads to the ground state of 176 Lu. The data from Oklo show clearly that 176 Lu is 
depleted in the reactor zones. But as concluded in ref. [9j, the degree of depletion will be a reliable indicator of the 
temperature only if an improved measurement of £> 9 (175) = a 9 /(a 9 + a m ) is performed, and if alternate explanations 
@ for depletion are ruled out. 

One alternate explanation for 176 Lu depletion lies in the possibility of processing of the lutetium isotopes in Oklo 
due to photo excitation of the isomeric state in 176m Lu by 176 Lu(7,7') fluorescence. Such a process is well known 
in astrophysics and is an important channel for burning 176 Lu in stellar environments [ill The isomeric state 
decays to 176 Hf with a half-life of 3.6 hr, and therefore provides an alternate path for removing 176 Lu. Here we explore 
whether this could have been be an effect in the 7-ray fluxes in the reactors, taking advantage of newly developed 
data libraries for fission decay chains. 

In 235 TJ thermal neutron fission, about 6.6 MeV is released in the form of prompt 7 rays, about 6.5 MeV as f3 rays, 
and about 6.3 MeV as 7 rays following f3 emission. The sum of the delayed (3 and 7 energy released during the decay 
of fission products is called decay heat and varies as a function of time, f(t), after a single fission event at t=0. 

Beginnin g w ith the work of Way and Wigner (l3| . many calculations and measurements of f(t) have been performed, 
see Tobias [14| and Dickens [l5j for reviews. Typically, f(t) ~ is found for times greater than several seconds 
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1 . With the development of more comprehensive nuclear data libraries based on level schemes derived from high 
resolution Ge-detector data, summation method calculations became widely accepted. These calculations gave good 
agreement with measurements except at the shorter times associated high Q-value f3 decays. High Q-values feed 
levels at high excitation that can decay by emission of weak and (or) high energy 7 rays easily missed in Ge-detector 
measurements. As a result, the libraries were incomplete. Recently, Total Absorption Gamma Spectrometer (TAGS) 
data have been included in libraries. This has eliminated the discrepancies. In particular, Algora et al. [13] were able 
to report 239 Pu decay heat calculations in excellent agreement with experiments for shorter times. 



We follow this approach for calculating 235 U decay heat 7-spectra |18|, converting to fluxes using standard 
energy deposition conversion coefficients. We apply Monte Carlo modeling to calculate the prompt 7-ray flux in 
the Oklo reactor zone. With these fluxes in hand, and a model of how the reactor operated, we can then estimate 
photo-excitation constants A 7i7 / for burning 176 Lu through the isomeric state 176m Lu. 



2. Prompt fission 7-ray flux in Oklo reactor zone RZ10. 

The MCNP code allows modeling of neutron transport and also provides the energy dependence of the prompt 
7-ray flux. We use the same input for the Oklo zone RZIO as in our previous work The model of a reactor zone 
is a flat cylinder of 70 cm height, 6 m diameter, surrounded by a I m thick reflector consisting of water saturated 
sandstone. As for any reactor, Oklo criticality is determined by the geometry and the composition of the active zone. 
Oklo reactor zones include uraninite UO2, gangue (oxides of different metals with water of crystallization) and water. 
The total density of the active core material at ancient times was about 3 g cm~ 3 for RZIO with only 30 wt. % of 
UO2 in the RZIO dry ore. The hydrogen to uranium atomic ratio in our model was jf 1 - — 13.0 and the multiplication 
coefficient of the fresh core was k e j^=1.036. Detailed composition and neutronic parameters of the RZIO reactor zone 
are given in [j| . The rate of fission can be deduced from the Hidaka and Holliger model [2(3] , which found an average 
RZIO neutron fluence of 0.65 kb _1 over a time duration of 160 kyr. 

From analysis of xenon isotope abundances in Oklo grains of aluminum phosphate, Meshik et al. (2lj concluded 
the reactors operated cyclically, with reactor-on periods of about 0.5 hr (1800 s), separated by dormant reactor-off 
periods of 2.5 hr (9000 s). We use this periodicity in calculating absolute fluxes. Adapting to the pulse mode, we 
find 1.03 • 10 18 fissions for the 1800 s the reactor is on, with a thermal power of about 18 kW. We include only 235 U 
fission here since 238 U and 239 Pu fission was found by Hidaka and Holliger to contribute less than ten percent to the 
neutron yields. 

The prompt 7-ray flux during the reactor-on period is shown as the upper line in Fig. 1. The other flux shown 
refers to the delayed heat and is discussed in the next sections. The total prompt 7-ray flux is about 3 ■ 10 9 7 cm" 2 
s _1 . The electromagnetic energy per fission corresponds to about 14.7 MeV, greater than the 6.6 MeV of prompt 
7-ray fission energy due to neutron capture on the materials of the reactor. 

3. Decay heat 7-ray spectra from fission 

A large number of radioactive nuclides are produced from the fission of an actinide target such as 235 U. Of relevance 
in describing the time evolution of the decay network are the longer-lived levels, ground states and isomers, here 
called materials. While the reactor is operating, the materials population satisfies the following set of linearly-coupled 
differential equations: 



dN t {t) 
dt 



-XiNi(t) + J2 *kPkiN k (t) + rFYi, (1) 



where the decay constant of the material i is Ai = ^n(2)/t 1 / 2 j, ti/2i is the half-life for the material i, Phi is the 
probability that the material k will populate in its decay the material i, r is the fission rate and FYi is the independent 
fission yield for the material i. When the reactor is not operating, the equations for iVj(T) are the same but with 
r=0. 

For each material, we have obtained the Ge-detector 7-ray spectrum Ii(E 7 ) from the ENDF/B-VII.l decay data sub- 
library [23|], using a detector resolution of 2 keV (FWHM) for discrete lines and a A_E 7 =0.5 keV binning. The spectrum 
is defined so that Ii(i? 7 )A£ , 7 gives the absolute probability of observing 7-rays with energies in the (i? 7 — 0.5Ai? 7 , 
E 1 + 0.5A.E 7 ) interval per decay of the material i. 

During the period t, t + St, the 7-ray spectrum obtained by adding the contribution of all the radioactive nuclides 
in the core is given by: 



1 Note this is not the same as the time dependence of decay heat following shutdown of a long-running reactor [l7j| . 
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I(£ 7 , t, t + St) = J2 XiNi(t)Ii(E^5t, (2) 

i 

which can be integrated over time to obtain: 

I(K 1 ,t ,h)=y2li(E 1 ) XiNi(t)dt. (3) 

i J to 

The mean electromagnetic energy (EEM), and average 7 energy of the spectrum are calculated as: 



EEM(t ,t 1 ) = J I(E 7 ,t ,h)E 7 dE^, (4) 

< E 7 (t ,h) >= J I{E 1 ,t ,ti)E 1 dE 1 / J I(E 1 ,t ,t 1 )dE 1 . (5) 

Equations (Q]) were solved numerically using the fission yields from the JEFF-3.1 library [22j and the decay data 
from the ENDF/B- VII. 1 library [13. 

To confirm the correctness of our procedure, we first calculated delayed 7-ray spectra for time intervals £=1-1800 
s and i=l-9000 s following a single fission event at t~0. These time intervals match the cycling times in our reactor 
model. Our EEM values are shown in the second column of Table 1. In columns three through five, we compare to 
values obtained by numerical integration of f(t) as given in Refs [HI H2, Hlj] ■ The last column is an average of these 
three reference values. Our EEM values are about 0.40 MeV lower because they are derived from point-wise 7-ray 
spectra. If we use the EEM values in the ENDF/B-VII.l library, which includes the latest TAGS measurements, 
we obtain EEM (1, 1800)=4.03 MeV and EEM(1, 9000)=4.77 MeV, in agreement with the average of other results. 
Because of the pandemonium effect j24|, it is a well-known fact that one obtains a lower EEM value when using 
point-wise spectra, in this case about 10%. 



TABLE I: Electromagnetic energies EEM(T) (in MeV per fission) in decay heat of 235 U. 



Time interval, s 


Present work 


ENDF/B-VII.F23] 


Dickensfl5] 


JEFF-3.1 [22] 


Average value 


1-1800 


3.60 


4.12 


4.10 


4.00 


4.07 


1-9000 


4.34 


4.84 


4.74 


4.70 


4.76 



We are interested in the 7-ray flux in the pulsed cycling mode of the reactor operation, not just for a fission event at 
t = 0. We therefore calculated next the 7-ray spectra during the 0.5 hr pulse (1800 s) and during the 2.5 hr (9000 s) 
cooling period, assuming one fission per second during the reactor-on pulse. In these calculations we took A^=0, that 
is, a fresh core. The spectra for a 0.15 to 10 MeV energy range are shown in Figs. 2 and 3, the former on a log scale, 
and the latter on a linear scale to show the fine structure of the lower energy portion of the spectrum. The EEM 
value for the reactor-on 0-1800 s spectrum is 3.49 MeV/fission, with < E 1 > equal to 0.732 MeV. For the reactor-off 
1800-10800 s spectrum, EEM is equal to 1.24 MeV/fission and < E 1 > is slightly higher, 0.813 MeV. Adding the 
spectra gives EEM(0, 10800) = 4.73 MeV. 

Using these EEM values and knowing the number of fissions per pulse, we can now calculate the prompt, delayed 
and equilibrium components of the electromagnetic energy for the reactor-on and reactor-off time intervals. These 
values are shown in columns two through four of Table II. The prompt entry corresponds to the energy from prompt 
fission and neutron capture 7-rays. This is zero when the reactor is off. The delayed entries are derived from the 7-ray 
decay EEM values given earlier. Both prompt and delayed components originate from one single reactor-on pulse 
of power 18 kW. In addition to these components, there will also be equilibrium decay heat associated with fission 
from the N previous pulses (N >>1). For our purposes this can be estimated over short periods of time simply by 
using the standard relation for the power of the after heat compared to the power of the reactor: P/ lea t = 0.066P atI 
[ltij . Taking the average thermal power in our model to be Pan =3 kW, and noting only half the energy in the after 
heat is electromagnetic, we then have P 7 = 0.033P atI , which leads to the values shown in column four. We see the 
equilibrium flux within the pulse will be a factor of 1.12/3.59 = 0.31 smaller than the delayed flux within the pulse. 

The last column is the energy release per unit time summed over all electromagnetic components. Assuming similar 
spectral shapes of the components, these latter quantities will be proportional to the total 7-ray flux during the 
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TABLE II: Electromagnetic energies (in 10 18 MeV) for the pulsed model of reactor RZ10. 



AT, s 


prompt 


delayed 


equilibrium 


energy rate, 10 ib MeV/s 


0-1800 


15.14 


3.59 


1.12 


11.03 


1800-10800 




1.28 


5.62 


0.77 



reactor-on and reactor-off pulses. We conclude that the flux in time intervals between pulses (1800—10800 s) will be 

about 7% of the flux during the 1800-s reactor-on pulse. 

4. Conversion of the delayed 7-ray spectra to 7-ray fluxes 

To estimate the photo-excitation parameters A 7i7 < in 176 Lu it is necessary to know 7-ray fluxes, not simply Ge- 
detector spectra. The prompt flux is obtained directly as output of MCNP 19]. We get the delayed flux using the 
7-ray dose to photon fluence conversion factors published by the International Commission on Radiological Protection 
(ICRP) [25j]. Our conversion assumes that all 7-rays produced in Oklo reactors are absorbed by the active core 
materials, both fertile and nonfertile. 

The conversion coefficients k(E), are listed in Table III, where the absorbed dose is in units of erg g _1 and the 
photon fluence is in units of 7 cm~ 2 . A polynomial fit to these data gives k(E) = 0.114_E 3 — 1.091.E 2 + 5.618-E — 0.189 
with E in MeV. 



TABLE III: Fluence-to-dose conversion coefficients k(E-y) for 7-rays with energy _B 7 [25[. 



E-,, keV 


100 


200 


300 


400 


500 


600 


800 


1000 


1500 


2000 


3000 


4000 


k(Ey), 10" B erg g-77 cm" 2 


0.37 


0.86 


1.38 


1.89 


2.38 


2.84 


3.69 


4.47 


6.14 


7.55 


9.96 


12.10 



For the Oklo zone RZ10, the total mass is 60 tonne. To convert the reactor-on spectrum of Fig. 2 to a flux, we 
therefore multiply by the energy E~, divide by k(Ej), divide by the total mass of the reactor zone and the 1800 
s accumulation time, and normalize to the 235 U fission rate 5.71 • 10 14 s _1 . Wc multiply this flux by 1.31 to take 
into account the contribution of the equilibrium flux. The resulting reactor-on delayed + equilibrium flux, binned in 
100-keV intervals, is shown in Fig 1 and can now be compared to the prompt flux. We see it is typically an order of 
magnitude smaller. However, at a few energies it does actually exceed the prompt flux. 

5. Implications of 7-ray flux estimates for 176 Lu/ 175 Lu thermometry. 

Lutetium thermometry is based on the dependence of the 176 Lu/ 175 Lu isotopic ratio on the operating temperature of 
a reactor in which 176 Lu is burned (and partially restituted) by neutron capture reactions. As detailed in Ref. 0, this 
process is not described by a single exponential decay constant. However, to set a time scale for judging the impact 
of the 7-ray flux, we can introduce an effective constant AJj" = 3.7 • 10~ 13 s _1 (effective half-life of 60 kY) based on 
the factor 6.4 reduction in the isotopic ratio over the 160 kY operating time of RZ10 [|. Then we write the total 
constant for disappearance of 176 Lu as A = A„ + A 7 . 7 ', where the second term corresponds to the photo-excitation 
process as an alternate explanation for 176 Lu depletion. We omit the beta decay constant for 176 Lu Ap = 5.8 • 10 -19 
s _1 because it is by six order of magnitude less than A^. 

Experimental data [ill . IT2I ] confirm that long lived 176 Lu in the photon bath of celestial bodies can be partially 
transformed into metastable 176m Lu by photons with energies around 880, 1060, 1330, and 1660 keV. Higher energy 
photons may also contribute [26| . These photon energies correspond to excited states of 176 Lu with specific spins and 
parities which allow them to act as mediators for photo excitation of the isomeric state. The rate A 7 . 7 '(i?i) of photo 
excitation of the isomeric state 176 Tu in photon inelastic scattering through an intermediate state (IS) at energy Ei 
is given by 

A 7) y(£i)= f ^{E)a ltY (E,E i )dE=^{E i )a™ t 1 ,(E i ). (6) 

Here $ 7 (i?i) is the differential photon flux at energy Ei having units of keV -1 cm~ 2 s -1 , and o'^iiEi) — 
J er 7 . 7 <(-E, Ei)dE is the integrated cross section for the IS. When several IS contribute, A 7i7 ' will be a sum over 
the individual IS contributions. 

In writing the rate on the right as the product of two factors we are implicitly assuming the flux is continuous, 
and varying slowly over the narrow resonance energy Ei. This is true in stellar environments, and for example, in 
bremsstrahlung experiments. It is not necessarily the case in our situation where the photon spectra are made up from 
a sum over many discrete 7-ray lines. However, our fluxes propagate in the dense environment of the Oklo reactors, 
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and as discussed by von Neumann-Cosel et al. [27| . Compton scattering can broaden otherwise discrete spectra quite 
significantly. The MCNP code models Compton scattering, pair-production and electron bremsstrahlung fully, leading 
to the continuous prompt spectrum shown in Fig. 1 . The delayed spectra show more structure, but for purposes of 
estimating upper bounds on the photo-excitation process we assume our 100 keV averaging procedure will serve as a 
useful approximation. 

Reported values of the photo excitation cross sections given in the literature include (-E 7 )=33.4 mb-eV for 
i? 7 =839 keV (this is an upper limit in Ref. [12j . for a laboratory environment, as opposed to a fully ionized envi- 
ronment), and higher values tr™y (E 7 )=140 mb-keV and 350 mb-keV, for 4- and 6-MeV bremsstrahlung irradiations 
respectively, with an assumed IS energy of 2.125 MeV [26j . 

With these cross sections, and the fluxes $ 7 (-E) of Fig. 1 at hand, we can now calculate A 7j7 < for Oklo reactor 
RZ10. Taking into account the pulse structure, we use a weighted average flux consisting of 1/6 of the sum of the 
prompt, delayed and equilibrium fluxes while the reactor is on, and 5/6 of 7% of this sum while the reactor is off. 

The option of £^=839 keV, with a total spectral flux of 0.86 • 10 6 7 cm" 2 s _1 keV -1 , leads to A 7) y = 0.26 • 10" 22 
s _1 which is negligible. The 6-MeV bremsstrahlung option with E 1 = 2.1 MeV, has a total spectral flux of 1.9 ■ 10 5 7 
cm~ 2 s _1 keV" 1 and gives A 7i7 ' = 6.8 • 10 -20 s _1 . However, even this value is less than the \%ff by seven orders of 
magnitude. We conclude that destruction of 176 Lu in the Oklo reactors is not influenced by any part of the reactor-on 
photon flux. 



After complete shutdown, the power of the 7-ray decay heat decreases approximately as t~ 3 [Tq, reaching about 

1% of the equilibrium value after a year. At this point it can be ignored. We conclude therefore that the decay heat 
is also not able to change the lutetium isotopic ratio even over the long period the reactor has been shut down. 

The particular values we have found are specific to the uniform pulsed mode of operation that we have as- 
sumed. However, the burn rate for steady state reactor operation will not be much different since our result is 
determined mainly by the prompt flux, which scales inversely with the live time of the reactor. Absent a many 
orders of magnitude larger photo-excitation cross section through as yet undetermined levels, we see the photon 
intensities in the reactor are insufficient to alter the 176 Lu/ 175 Lu isotopic ratios associated with neutron transmutation. 

6. Conclusions 

We have developed realistic models of the prompt and delayed 7-ray fluxes in Oklo natural nuclear reactors, taking 
advantage of recent releases of the databases and methodology that accurately represent short time decay heat in 
fission processes. We have compared 176 Lu transmutation rates associated with photo excitation to transmutation 
rates associated with neutron capture processes. In contrast to astrophysical processes, we find 176 Lu/ 175 Lu isotopic 
ratios in Oklo are not influenced by decay heat electromagnetic radiation, either during reactor operation, or after 
reactor shutdown. Lutetium thermometry, as recently studied [9j , is therefore applicable to analyses of Oklo reactor 
data. 
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FIG. 1: (Color online). Prompt and delayed 7-ray fluxes $ 7 (i5) in Oklo reactor zone RZ10 during the period the reactor is 
on. The fluxes are calculated for an 18 kW reactor cycling on for 0.5 hr, and off for 2.5 hr. The prompt flux, from MCNP, 
is the upper line. The lower line is the delayed flux for one 0.5 hr fresh-core reactor-on pulse, multiplied by 1.31 to take into 
account the equilibrium flux associated with the N previous reactor-on pulses (N >>1). The statistical uncertainty in the 
prompt spectrum simulation is 5 %. The structure in the delayed flux is due to incomplete averaging of contributions from the 
discrete lines shown in Fig. 3. 



8 




2000 4000 6000 8000 10000 



Energy (keV) 



FIG. 2: (Color online). Time-integrated 7-ray spectra in the 0.15 MeV to 10 MeV range for the T = — 1800 s reactor-on 
pulse (upper) at one fission/sec, and for the reactor-off period T = 1800 — 10800 (lower), the latter scaled down by a factor of 
10 3 for clarity. Counts are per 1-eV energy interval. 
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FIG. 3: (Color online). Time-integrated 7-ray spectra in the 0.15 MeV to 5 MeV range for the T = — 1800 s reactor-on pulse 
(upper), and for the reactor-off period T = 1800 — 10800 (lower). Counts are per 1-eV energy interval. The discrete nature 
of the spectra is more evident as compared to Fig. 2, but (see text) can still be treated as quasi-continuous for purposes of 
estimating photo-excitation probabilities. 



